Development of single nucleotide polymorphism (SNP) markers from the mango (Mangifera indica) transcriptome for mapping and estimation of genetic diversity

The development of molecular markers for genomic studies in Mangifera indica (mango) will allow marker-assisted selection and identification of genetically diverse germplasm, greatly aiding mango breeding programs. We report here our identification of thousands of unambiguous molecular markers that can be easily assayed across genotypes of the species. With origin centered in Southeast Asia, mangos are grown throughout the tropics and subtropics as a nutritious fruit that exhibits remarkable intraspecific phenotypic diversity. With the goal of building a high density genetic map, we have undertaken discovery of sequence variation in expressed genes across a broad range of mango cultivars. A transcriptome sequence reference was built de novo from extensive sequencing and assembly of RNA from cultivar 'Tommy Atkins'. Single nucleotide polymorphisms (SNPs) in protein coding transcripts were determined from alignment of RNA reads from 24 mango cultivars of diverse origins: 'Amin Abrahimpur' (India), 'Aroemanis' (Indonesia), 'Burma' (Burma), 'CAC' (Hawaii), 'Duncan' (Florida), 'Edward' (Florida), 'Everbearing' (Florida), 'Gary' (Florida), 'Hodson' (Florida), 'Itamaraca' (Brazil), 'Jakarata' (Florida), 'Long' (Jamaica), 'M. Casturi Purple' (Borneo), 'Malindi' (Kenya), 'Mulgoba' (India), 'Neelum' (India), 'Peach' (unknown), 'Prieto' (Cuba), 'Sandersha' (India), 'Tete Nene' (Puerto Rico), 'Thai Everbearing' (Thailand), 'Toledo' (Cuba), 'Tommy Atkins' (Florida) and 'Turpentine' (West Indies). SNPs in a selected subset of protein coding transcripts are currently being converted into Fluidigm assays for genotyping of mapping populations and germplasm collections. Using an alternate approach, SNPs (144) discovered by sequencing of candidate genes in 'Kensington Pride' have already been converted and used for genotyping.

The use of high resolution melting (HRM) to map single nucleotide polymorphism markers linked to a covered smut resistance gene in barley

Using an established genetic map, a single gene conditioning covered smut resistance, Ruh.7H, was mapped to the telomere region of chromosome 7HS in an Alexis/Sloop doubled haploid barley population. The closest marker to Ruh.7H, abg704 was 7.5 cM away. Thirteen loci on the distal end of 7HS with potential to contain single nucleotide polymorphisms (SNPs) were identified by applying a comparative genomics approach using rice sequence data. Of these, one locus produced polymorphic co-dominant bands of different size while two further loci contained SNPs that were identified using the recently developed high resolution melting (HRM) technique. Two of these markers flanked Ruh.7H with the proximal marker located 3.8 cM and the distal marker 2.7 cM away. This is the first report on the application of the HRM technique to SNP detection and to rapid scoring of known cleaved amplified polymorphic sequence (CAPS) markers in plants. This simple, precise post-PCR technique should find widespread use in the fine-mapping of genetic regions of interest in complex cereal and other plant genomes.

Characterization of disease resistance gene candidates of the nucleotide binding site (NBS) type from banana and correlation of a transcriptional polymorphism with resistance to Fusarium oxysporum f.sp. cubense race 4

Most plant disease resistance (R) genes encode proteins with a nucleotide binding site and leucine-rich repeat structure (NBS-LRR). In this study, degenerate primers were used to amplify genomic NBS-type sequences from wild banana (Musa acuminata ssp. malaccensis) plants resistant to the fungal pathogen Fusarium oxysporum formae specialis (f. sp.) cubense (FOC) race 4. Five different classes of NBS-type sequences were identified and designated as resistance gene candidates (RGCs). The deduced amino acid sequences of the RGCs revealed the presence of motifs characteristic of the majority of known plant NBS-LRR resistance genes. Structural and phylogenetic analyses grouped the banana RGCs within the non-TIR (homology to Toll/interleukin-1 receptors) subclass of NBS sequences. Southern hybridization showed that each banana RGC is present in low copy number. The expression of the RGCs was assessed by RT-PCR in leaf and root tissues of plants resistant or susceptible to FOC race 4. RGC1, 3 and 5 showed a constitutive expression profile in both resistant and susceptible plants whereas no expression was detected for RGC4. Interestingly, RGC2 expression was found to be associated only to FOC race 4 resistant lines. This finding could assist in the identification of a FOC race 4 resistance gene.

The assessment of TGGE for the detection of interspecific and intergeneric DNA-marker polymorphism within Solanaceae

RFLP markers are currently the most appropriate marker system for the identification of uncharacterised polymorphism at the interspecific and intergeneric level. Given the benefits of a PCR-based marker system and the availability of sequence information for many Solanaceous cDNA clones, it is now possible to target conserved fragments, for primer development, that flank sequences possessing interspecific polymorphism. The potential outcome is the development of a suite of markers that amplify widely in Solanaceae. Temperature gradient gel electrophoresis (TGGE) is a relatively inexpensive gel-based system that is suitable for the detection of most single-base changes. TGGE can be used to screen for both known and unknown polymorphisms, and has been assessed here, for the development of PCR-based markers that are useful for the detection of interspecific variation within Solanaceae. Fifteen markers are presented where differences between Lycopersicon esculentum and L. pennellii have been detected by TGGE. The markers were assessed on a wider selection of plant species and found to be potentially useful for the identification of interspecific and intergeneric polymorphism in Solanaceous plants.

Finding genes for economically important traits: Brahman cattle puberty.

Age at puberty is an important component of reproductive performance in beef cattle production systems. Brahman cattle are typically late-pubertal relative to Bos taurus cattle and so it is of economic relevance to select for early age at puberty. To assist selection and elucidate the genes underlying puberty, we performed a genome-wide association study (GWAS) using the BovineSNP50 chip (similar to 54 000 polymorphisms) in Brahman bulls (n = 1105) and heifers (n = 843) and where the heifers were previously analysed in a different study. In a new attempt to generate unbiased estimates of single-nucleotide polymorphism (SNP) effects and proportion of variance explained by each SNP, the available data were halved on the basis of year and month of birth into a calibration and validation set. The traits that defined age at puberty were, in heifers, the age at which the first corpus luteum was detected (AGECL, h(2) = 0.56 +/- 0.11) and in bulls, the age at a scrotal circumference of 26 cm (AGE26, h(2) = 0.78 +/- 0.10). At puberty, heifers were on average older (751 +/- 142 days) than bulls (555 +/- 101 days), but AGECL and AGE26 were genetically correlated (r = 0.20 +/- 0.10). There were 134 SNPs associated with AGECL and 146 SNPs associated with AGE26 (P < 0.0001). From these SNPs, 32 (similar to 22%) were associated (P < 0.0001) with both traits. These top 32 SNPs were all located on Chromosome BTA 14, between 21.95 Mb and 28.4 Mb. These results suggest that the genes located in that region of BTA 14 play a role in pubertal development in Brahman cattle. There are many annotated genes underlying this region of BTA 14 and these are the subject of current research. Further, we identified a region on Chromosome X where markers were associated (P < 1.00E-8) with AGE26, but not with AGECL. Information about specific genes and markers add value to our understanding of puberty and potentially contribute to genomic selection. Therefore, identifying these genes contributing to genetic variation in AGECL and AGE26 can assist with the selection for early onset of puberty.

Association genetics in Corymbia citriodora subsp variegata identifies single nucleotide polymorphisms affecting wood growth and cellulosic pulp yield

Wood is an important biological resource which contributes to nutrient and hydrology cycles through ecosystems, and provides structural support at the plant level. Thousands of genes are involved in wood development, yet their effects on phenotype are not well understood. We have exploited the low genomic linkage disequilibrium (LD) and abundant phenotypic variation of forest trees to explore allelic diversity underlying wood traits in an association study. Candidate gene allelic diversity was modelled against quantitative variation to identify SNPs influencing wood properties, growth and disease resistance across three populations of Corymbia citriodora subsp. variegata, a forest tree of eastern Australia. Nine single nucleotide polymorphism (SNP) associations from six genes were identified in a discovery population (833 individuals). Associations were subsequently tested in two smaller populations (130160 individuals), validating our findings in three cases for actin 7 (ACT7) and COP1 interacting protein 7 (CIP7). The results imply a functional role for these genes in mediating wood chemical composition and growth, respectively. A flip in the effect of ACT7 on pulp yield between populations suggests gene by environment interactions are at play. Existing evidence of gene function lends strength to the observed associations, and in the case of CIP7 supports a role in cortical photosynthesis.

The use of genetic correlations to evaluate associations between SNP markers and quantitative traits

Open-pollinated progeny of Corymbia citriodora established in replicated field trials were assessed for stem diameter, wood density, and pulp yield prior to genotyping single nucleotide polymorphisms (SNP) and testing the significance of associations between markers and assessment traits. Multiple individuals within each family were genotyped and phenotyped, which facilitated a comparison of standard association testing methods and an alternative method developed to relate markers to additive genetic effects. Narrow-sense heritability estimates indicated there was significant additive genetic variance within this population for assessment traits ( h ˆ 2 =0.28to0.44 ) and genetic correlations between the three traits were negligible to moderate (r G = 0.08 to 0.50). The significance of association tests (p values) were compared for four different analyses based on two different approaches: (1) two software packages were used to fit standard univariate mixed models that include SNP-fixed effects, (2) bivariate and multivariate mixed models including each SNP as an additional selection trait were used. Within either the univariate or multivariate approach, correlations between the tests of significance approached +1; however, correspondence between the two approaches was less strong, although between-approach correlations remained significantly positive. Similar SNP markers would be selected using multivariate analyses and standard marker-trait association methods, where the former facilitates integration into the existing genetic analysis systems of applied breeding programs and may be used with either single markers or indices of markers created with genomic selection processes.

Differentiation of Indian, East Timorese, Papuan and Floridian Candidatus Liberibacter asiaticus' isolates on the basis of simple sequence repeat and single nucleotide polymorphism profiles at 25 loci

Japanese isolates of Candidatus Liberibacter asiaticus have been shown to be clearly differentiated by simple sequence repeat (SSR) profiles at four loci. In this study, 25 SSR loci, including these four loci, were selected from the whole-genome sequence and were used to differentiate non-Japanese samples of Ca. Liberibacter asiaticus (13 Indian, 3 East Timorese, 1 Papuan and 8 Floridian samples). Out of the 25 SSR loci, 13 were polymorphic. Dendrogram analysis using SSR loci showed that the clusters were mostly consistent with the geographical origins of the isolates. When single nucleotide polymorphisms (SNPs) were searched around these 25 loci, only the upstream region of locus 091 exhibited polymorphism. Phylogenetic tree analysis of the SNPs in the upstream region of locus 091 showed that Floridian samples were clustered into one group as shown by dendrogram analysis using SSR loci. The differences in nucleotide sequences were not associated with differences in the citrus hosts (lime, mandarin, lemon and sour orange) from which the isolates were originally derived.

High resolution mapping of Dense spike-ar (dsp.ar) to the genetic centromere of barley chromosome 7H

Spike density in barley is under the control of several major genes, as documented previously by genetic analysis of a number of morphological mutants. One such class of mutants affects the rachis internode length leading to dense or compact spikes and the underlying genes were designated dense spike (dsp). We previously delimited two introgressed genomic segments on chromosome 3H (21 SNP loci, 35.5 cM) and 7H (17 SNP loci, 20.34 cM) in BW265, a BC7F3 nearly isogenic line (NIL) of cv. Bowman as potentially containing the dense spike mutant locus dsp.ar, by genotyping 1,536 single nucleotide polymorphism (SNP) markers in both BW265 and its recurrent parent. Here, the gene was allocated by high-resolution bi-parental mapping to a 0.37 cM interval between markers SC57808 (Hv_SPL14)-CAPSK06413 residing on the short and long arm at the genetic centromere of chromosome 7H, respectively. This region putatively contains more than 800 genes as deduced by comparison with the collinear regions of barley, rice, sorghum and Brachypodium, Classical map-based isolation of the gene dsp.ar thus will be complicated due to the infavorable relationship of genetic to physical distances at the target locus.

Mapping the carriage of flaA-restriction fragment length polymorphism Campylobacter genotypes on poultry carcasses through the processing chain and comparison to clinical isolates

Poultry are considered a major source for campylobacteriosis in humans. A total of 1866 Campylobacter spp. isolates collected through the poultry processing chain were typed using flaA-restriction fragment length polymorphism to measure the impact of processing on the genotypes present. Temporally related human clinical isolates (n = 497) were also typed. Isolates were obtained from whole chicken carcass rinses of chickens collected before scalding, after scalding, before immersion chilling, after immersion chilling and after packaging as well as from individual caecal samples. A total of 32 genotypes comprising at least four isolates each were recognised. Simpson's Index of Diversity (D) was calculated for each sampling site within each flock, for each flock as a whole and for the clinical isolates. From caecal collection to after packaging samples the D value did not change in two flocks, decreased in one flock and increased in the fourth flock. Dominant genotypes occurred in each flock but their constitutive percentages changed through processing. There were 23 overlapping genotypes between clinical and chicken isolates. The diversity of Campylobacter is flock dependant and may alter through processing. This study confirms that poultry are a source of campylobacteriosis in the Australian population although other sources may contribute.

Genome-wide association analysis and differential expression analysis of resistance to Sclerotinia stem rot in Brassica napus

Brassica napus is one of the most important oil crops in the world, and stem rot caused by the fungus Sclerotinia sclerotiorum results in major losses in yield and quality. To elucidate resistance genes and pathogenesis-related genes, genome-wide association analysis of 347 accessions was performed using the Illumina 60K Brassica SNP (single nucleotide polymorphism) array. In addition, the detached stem inoculation assay was used to select five highly resistant (R) and susceptible (S) B. napus lines, 48 h postinoculation with S. sclerotiorum for transcriptome sequencing. We identified 17 significant associations for stem resistance on chromosomes A8 and C6, five of which were on A8 and 12 on C6. The SNPs identified on A8 were located in a 409-kb haplotype block, and those on C6 were consistent with previous QTL mapping efforts. Transcriptome analysis suggested that S. sclerotiorum infection activates the immune system, sulphur metabolism, especially glutathione (GSH) and glucosinolates in both R and S genotypes. Genes found to be specific to the R genotype related to the jasmonic acid pathway, lignin biosynthesis, defence response, signal transduction and encoding transcription factors. Twenty-four genes were identified in both the SNP-trait association and transcriptome sequencing analyses, including a tau class glutathione S-transferase (GSTU) gene cluster. This study provides useful insight into the molecular mechanisms underlying the plant's response to S. sclerotiorum.

Identification of I-7 expands the repertoire of genes for resistance to Fusarium wilt in tomato to three resistance gene classes

The tomato I-3 and I-7 genes confer resistance to Fusarium oxysporum f. sp. lycopersici (Fol) race 3 and were introgressed into the cultivated tomato, Solanum lycopersicum, from the wild relative Solanum pennellii. I-3 has been identified previously on chromosome 7 and encodes an S-receptor-like kinase, but little is known about I-7. Molecular markers have been developed for the marker-assisted breeding of I-3, but none are available for I-7. We used an RNA-seq and single nucleotide polymorphism (SNP) analysis approach to map I-7 to a small introgression of S. pennellii DNA (c. 210 kb) on chromosome 8, and identified I-7 as a gene encoding a leucine-rich repeat receptor-like protein (LRR-RLP), thereby expanding the repertoire of resistance protein classes conferring resistance to Fol. Using an eds1 mutant of tomato, we showed that I-7, like many other LRR-RLPs conferring pathogen resistance in tomato, is EDS1 (Enhanced Disease Susceptibility 1) dependent. Using transgenic tomato plants carrying only the I-7 gene for Fol resistance, we found that I-7 also confers resistance to Fol races 1 and 2. Given that Fol race 1 carries Avr1, resistance to Fol race 1 indicates that I-7-mediated resistance, unlike I-2- or I-3-mediated resistance, is not suppressed by Avr1. This suggests that Avr1 is not a general suppressor of Fol resistance in tomato, leading us to hypothesize that Avr1 may be acting against an EDS1-independent pathway for resistance activation. The identification of I-7 has allowed us to develop molecular markers for marker-assisted breeding of both genes currently known to confer Fol race 3 resistance (I-3 and I-7). Given that I-7-mediated resistance is not suppressed by Avr1, I-7 may be a useful addition to I-3 in the tomato breeder's toolbox.

SorGSD: a sorghum genome SNP database

Sorghum (Sorghum bicolor) is one of the most important cereal crops globally and a potential energy plant for biofuel production. In order to explore genetic gain for a range of important quantitative traits, such as drought and heat tolerance, grain yield, stem sugar accumulation, and biomass production, via the use of molecular breeding and genomic selection strategies, knowledge of the available genetic variation and the underlying sequence polymorphisms, is required.